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The present contribution deals with surface terms appearing immediately in distributions and 
partition functions of statistical ensembles. It is shown that all ensembles under study, including 
ordinary canonical and grand canonical ensembles, involve surface terms. For a canonical ensemble 
both surface and volume terms correspond to a closed system. For a grand canonical ensemble the 
volume term corresponds to an open system, while the surface one - to a closed one. Finally, for 
the recently introduced open statistical ensemble the specific feature of which is the consideration 
of some surrounding region both volume and surface terms correspond to an open system. 

In conclusion, surface particles at solid/fluid boundary are interpreted as particles corresponding 
to number density oscillations near the surface; this completely agrees with the earlier introduced 
concept of surface tension at such a boundary. 
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I. INTRODUCTION 

In paper [l[ an open statistical ensemble (OSE) is stud- 
ied, and comparison with a grand canonical ensemble 
(GCE) is made. The main peculiarity of OSE is correct 
consideration of surface terms. 

OSE approach enables one to examine surface terms 
directly on the basis of fundamental statistical relations 
rather than to introduce them additionally. 

Unlike GCE, for OSE the analogs of compressibility 
equation for higher order and 'mixed' correlations hold 
rigorously. 

The present paper shows that surface terms can be 
found even in a standard canonical ensemble (CE), and 
also correspond to the boundary between fluid and hard 
solid. 

Thermodynamic analysis verifies the form of surface 
terms in CE and the expression for fi-potential of the 
CE as a closed system with the surface. 

In papers i and GCE surface terms are examined. 
Comparison between OSE and GCE made in l| demon- 
strates that surface terms are equal in value and opposite 
in sign. In this contribution we reveal the meaning of this 
fact related to hybrid nature of GCE. 

Paper [l[ shows that unlike equilibrium two-phase sys- 
tem, for the case at hand the coefficient of surface tension 
depends on two variables, namely, pressure and tempera- 
ture, and corresponds to the boundary between fluid and 
hard solid. This results in the non-zero number of sur- 
face particles for such a boundary. This paper establishes 
their physical meaning as particles lying within the limits 
of near-surface oscillations or other near-surface distor- 
tions corresponding to local deviation of density from the 
mean one. 



II. PRIMARY DEFINITIONS 
A. Canonical ensemble 

Let us start with configuration integral 

Z N = J exp(-@u£ N )dri...dr N . 



(1) 



Here N is the number of particles in the system, /3 = 
1/ksT, ks is the Boltzmann constant, T is the temper- 
ature, N is the interaction energy of particles, V is 
the system volume. Integration is performed over the 
coordinates of all particles of the ensemble ri...r^. 

The probability density of finding a given space config- 
uration of a particular set of particles [3|, p. 181] is given 
by the expression 

P[ k \ = ir\ exp(-pul, N )dr k+1 ...dr N . (2) 

J 

V 

Passing to distribution functions for arbitrary set of par- 
ticles, we have 



0) _ 

Pc,l...k ~ 



(N - k) 



(fc) 



(3) 



where k defines the probability density of finding a 
given configuration of k arbitrary particles for the canon- 
ical ensemble. 
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B. Grand canonical ensemble 

Average equality (J3j over fluctuations of the number 
of particles, i.e., apply the operation Y^n=o -Pjv to both 
its sides. Here 



(4) 
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is the probability for GCE to have a certain number N 
of particles within the volume V. z is the activity 



(5) 



where fi is the chemical potential, A = h/\/27rmkBT, h is 
the Planck constant, m is the particle mass, and Sy is a 
large partition function of the system of the volume V 
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We obtain 
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J exp(-(3U^^ +k )dr k+1 ...dr k+N . 
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q g \ u specify the probability density of finding certain 
configuration of arbitrary particles for GCE. For ideal gas 
Qqi k = Q k , where g = N /V is the density of particles. 



others. We give several first functions for a homogeneous 
medium. 



(10) 
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Obviously, these functions can be formed for any ensem- 
ble, and the equality valid for arbitrary k > 1 may serve 
as their definition 

{n} Q— 1 



(11) 



Q = l 



where {«,} denotes some partitioning of a given set of 
k particles with the coordinates ri,...rk into disjoint 
groups {n a }, I is the number of groups of a specific par- 
titioning, k a is the size of the group with the number a, 
summation is taken over all possible partitionings. 
The inverse relation is valid 



C. Open statistical ensemble 



OSE functions corresponding to Qqi k are defined as 



(9) 



on condition that the region containing the coordinates 
of particles remains finite, and intensive parameters of 
the medium are kept unchanged. In practice, this means 
that we can use the definition of GCE correlation func- 
tions rather far away from the boundaries of the ensemble 
volume. 

In contrast to GCE correlation functions, the expan- 
sion Qi into a series in powers of activity involves no 
window functions, and integration in the integrals of the 
expansion coefficients is performed over the infinite space 



D. Localized correlation functions 

Application of logarithmic generating function gives 
rise to Ursell functions or localized correlation functions 
1,1- 

These functions have the property of locality: they 
decay rapidly when any groups of particles move apart 
including the case with one particle moving away from 
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with notation the same as in (fTTj) [6|, Chapter 9] . 



E. Ursell factors 

Ursell factors ti[ k \ are also called cluster functions. 
These are the functions that appear in the familiar pres- 
sure expansion in powers of activity (|3"T)) [1, p. 129], @, 
p.232]. 

They may be defined by the equality 



wffik = E^ 1 )'" 1 ^ - ^ II cxp(-/?£/ fc °({n Q })), 

{n} a=l 
I 

1 < k a <k, fc " = k, expi-fiU 1 ) = 1, (13) 

a=l 

where designations are the same as in (jlip. and the mean- 
ing of the condition exp(— /3J7 1 ) = 1 is evident. 
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As an example we give several first U\ k 

U[ 1] =1 (14) 
Ufl = exp(-pUl 2 ) - 1 

"$3 = ex P (-/3C/ 1 3 2i3 ) - eM-PUl 2 ) 

- exp(-/3f/ 1 2 3 ) - exp(-/3C/ 2 2 3 ) + 2 

Ursell factors also decay rapidly when any group of par- 
ticles, including a single particle, moves away. 

III. CANONICAL ENSEMBLE 

First we shall describe the mechanism of calculating 
the first surface terms of the series in activity in this 
ensemble, then general proof will be given. 

A. The series in the group size 

Let us employ expressions (fTTj) and (TT3"]) that yield ana- 
log (|12l) - the Bolzmann factor expansion in Ursell factors 

i 

exp(-/?C/ 1 fe .. fe )=^n w(fcQ) ({^})' 

{n} a=l 

I 

l<k a <k, E = k - ( 15 ) 

a=l 

As before, designations are the same as in (fTT|) . 



As far as we know, equation (fT5l) was first derived in 
4] , and was used for the terms of configuration integral 
logarithm expansion in powers of activity. 

So far it has probably remained unrealized that ex- 
pansion (|15[) may be treated as a series in the parameter 
which is the total number of particles in nonsingle groups 
(the sum of sizes of the partition groups except single 
particle groups). In view of 

UP = 1, (16) 

and changing the indices, we rewrite equation (| L5[) as 

JV I 

expR?^) = 1 + EE IT " ( * a) ({»«})> 

k=2 {n} a=l 
I 

2<k a <k, Y, k <* = k, (17) 

with notation also the same as in (|11[) except that parti- 
tionings are made over nonsingle groups. 

It is easily seen that, unlike (TI"5)) . (|T7|) explicitly demon- 
strates the absense of single particle groups (isolated par- 
ticles) . 

Thus the principal idea of this section is to apply (jTTJ) 
immediately to the entire Boltzmann factor included in 
Zjv, i.e., to use the partitioning of the whole combination 
of N particles contained in the closed system into groups. 

If expanded, (fT7|) has the form of a series in the total 
size of the groups 



exp(-^) = l+ £ U% + £ U% k + E Kk«+wS , wS+^MJ+Wi?w8} + .-. 
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Substitute p7|) in (fTJ). Making integration over the coor- 
dinates of free particles and grouping identical terms, we 
pass to a series for configuration integral 
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jrrij = k. 



(19) 



Here in the internal sum we go from the sum over parti- 
tionings to the sum over the partitioning topology: {m} 
is some partitioning of the group of k particles into sub- 



groups of the size j irrespective of the numbers of parti- 
cles, rrij is the quantity of subgroups of the size j, sum- 
mation is taken over all possible partitionings, and Cj has 
the form 



X...dTj. 



(20) 



Expression (|19p is similar to the expression given in 0, 
p. 128], but with single particle groups explicitly elimi- 
nated from it, and with traditional expression Vbj re- 
placed by Cj the meaning of which will become clear later 
on. In the expanded form 
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(3) 



1 



dr\...drA + . 



Generally speaking, this series (in density) is not a series 
in a small parameter, however, it is structurally similar 
to exponent expansion, so we can take a logarithm of it 
and see that the obtained series corresponds to correct 
volume and surface terms. In other words, when tak- 
ing a logarithm , for example, of the term with k = 2, 
we actually make allowance for the products of all pair 
groups: double, triple, etc., entering in sum (|2"Tj) . The 
analysis presented solves the problems with the series for 
configuration integral @, p. 227]. 



B. Separating out surface terms 

Introduce window functions of two types 

i> v {n) = $ ■■ 

and 
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(22) 



(23) 



Application of ipY and xY implies that integration is per- 
formed over the infinite space, unless otherwise specified. 

Rewrite expression (f2Tfl) with the use of these functions. 
This gives 



1 



k=l 



U^ j dr 1 ...drj. 



(24) 



The procedure of separating out surface terms [T| yields 

Cj = Vbj + Aaj, (25) 

where A is the surface surrounding the system, Vbj cor- 
responds to traditional definition [3j, p. 129] 



(26) 



and 



Aaj 



1 



J 



V 



i 



_fe=2 

or, with tj}Y replaced by 1 — \X between the brackets, 



U[ j) 1 dr 1 ...dr 



(27) 



Aaj = — 



xl) 



Lfe=2 



U[ i) 1 dr 1 ...dr r (28) 



After the product expansion the integrands in (|28p always 
involve at least one pair of particles on different sides of 
the system boundary, thus providing the proportionality 
of this expression to its surface pj. 

Below we show that the presence of surface terms in 
(|25p is not accidental, and that boundary mathematical 
effects coincide closely with physical ones. 

In pure form for a,j and bj we have 



U[ j) 1 dr 2 ...dr ] , 



(29) 



°i = J I i>i 



II(i-xD-i 



,fc=2 



Ui ' jdxidr2—drj, 



(30) 

where X\ is the coordinate perpendicular to the system 
surface, and the axis direction corresponds to positive 
increment. 

So the partition function of the canonical ensemble in- 
volves surface terms. Their exact form may be obtained 
by substituting (|23|) in (flT))) , and, if necessary, taking a 
logarithm of the series. 

The first terms of the expansion are 



Zn 



N\ (Vb 2 + Aa 2 ) 



(iV-2)! 



V 2 



N\ (Vb z + Aa 3 ) 
(N - 3)! V 3 



Nl [(Vbi + Aai) + \{Vb 2 +Aa 2 ) 



(N-4)\ 



V* 



(31) 



Taking a logarithm and in view of where Fn is free energy and £In is omega potential of 

Q the closed system, we get 

F N = Nk B T\n(NA 3 ) - Nk B T - k B T\nZ N ; 

n N = F N - fj,N, (32) n N = ~Vk B T(z + z 2 b 2 + z 3 b 3 + ...) 

-Ak B T{z 2 a 2 + z 3 a 3 + ...). (33) 
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Here we go from the series in powers of density to the 
series in powers of activity using the known expansion 



q(z) = z + z^2(n + l)z n b ri 

n=l 



(34) 



In the next section we shall establish the general form of 
the O-potential of the canonical ensemble, however, first 
it is necessary to consider grand canonical ensemble. 



IV. GRAND CANONICAL ENSEMBLE 

In [l[ the procedure of separating out volume and sur- 
face terms similar to that described above was applied to 
GCE partition function 



00 r 
t=i J 
It has been shown that 



IK 



U^} t dri...dr t . 



Z v = expp[VP(z,T)-A<T(z ) T)} ) 



(35) 



(36) 



where 



00 k r 

P{z,T) = zk B T + k B TY, 7 y U ( i k) . k dr 2 ...dr k , (37) 

fc=2 ' ^ 

is the well-known expansion for pressure, and for surface 
tension coefficient we have the expression 



a(z,T) = k B T^2 7f 



(38) 



x < 



i-f[(i-xJ 



U-f* t dxidr 2 --.dr t 



Equation (|36[) defines usual relation between partition 
function and GCE ^-potential 



- fc s TlnSy = fl = -PV + a A, 



(39) 



but the fi-potential itself has an unusual form. First, it 
corresponds to the system with a surface, second, now 
the surface tension coefficient depends on two variables. 
However, the first point is not surprising, since GCE is 
formed as a set of closed systems. Naturally, if we con- 
sider the systems each of which has a surface, this also 
refers to a combination of these systems. 

The surface term sign corresponds exactly to a closed 
system, and contradicts thermodynamic expression for 
fluctuation probability in an open one 

So GCE is an open system from the standpoint of vol- 
ume properties, since it describes adequately the fluctu- 
ation of the total number of particles. 

However, from the point of view of surface properties, 
it is a closed system. This is primarily supported by the 
fact that for GCE the mean number of particles is equal 



to the sum of volume and surface ones [l[. Probably, it 
is not coincidental that the author of paper 0] has come 
to the conclusion that GCE describes a "drop". 

Return to expression (flU)) . As mentioned above, it is 
actually identical to the expression derived in [3J, p. 128] 
with the substitution 



Vbj->Cj. 



(40) 



With saddle-point technique it has been proved [!, p. 150] 
that for the canonical ensemble we have the following ex- 
pression for the ri-potential accurate to the terms loga- 
rithmic in the number of particles 



(41) 



Proceeding as in [3j, p. 150] with the use of substitution 
(40), we obtain 



N 



-VP(z,T) + Aa(z,T), (42) 



j>i 



where P(z, T) and cr(z, T) are defined by expressions (|37p 
and (|38l) . respectively. 

Thus series (pl3"|) is nothing else but the first terms of 
the expansion in powers of activity of expression (|42l) . 



V. OPEN STATISTICAL ENSEMBLE 

First we give some principal results derived in [lj for 
OSE. Unlike ([35]), (|36| . for the partition function T v of 
this ensemble, we have the expressions 



*=i 



n 



Ui] t dr\...drt 



and 



T y = exp f3{VP(z,T) + Aa(z,T)}. 



(43) 



(44) 



In contrast to GCE, in OSE the mean number of particles 
contains no surface terms, thus the f2-potential of OSE - 
fiy - corresponds to a true open system 



Naturally, 



k B T\x\T\ 



-PV. 
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oA. 



(45) 



(46) 



Equality (|46[) seems logical, since the system thermody- 
namics is not to be determined solely by the probability 
of large deviations which is related to p$, and, therefore, 
to Ty- In other words, in (|4T)|) we see that, unlike (|39|) . 
the partition function of the system is determined both 
by thermodynamic potential (VL) related to the mean val- 
ues, and by the surface term related to large deviations. 
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VI. SURFACE PARTICLES 



So, since, according to (|38|) , a depends on two vari- 
ables, we obtain the non-zero value of the number of sur- 
face particles 



-A 



da 
Oft 



-0Az 



da 

Or: 



(47) 



or, using 



4=2 



i=2 



bti' t dri...dr t . 



(48) 

Equality (|4"5)l is almost apparent if we consider the famil- 
iar formula for density expansion into a series in powers 
of activity for GCE 



71=1 ' J 

and similar value for OSE 



'n+1 
.i=2 



Ui n . + nlidr 2 --dr n+1 . 

(49) 



g(z) = z + *Y, Z — U[ n+ n %dr 2 ...dr n+1 , (50) 

n=l n ' •* 

corresponding to (|57) . 



The number of surface particles in this case is identified 
with the value 



Ns = J [Qa(r,z) - g(z)} dr, 
v 



(51) 



and surface particle density gs — Ns /A is given by beau- 
tiful formula 



Lis 



[qg{x, z) - g{z)} dx, 



(52) 



it 



where before - the coordinate perpendicular to the 
surface, and L t - the transition region near the surface. 

In (fB"Tj) . (p)2")) we employ the expression for boundary 
density of particles at the outer boundary of GCE. How- 
ever, it is easily seen that the conclusion remains un- 
changed for the case of hard solid placed in the system. 
Moreover, result is the same for similar consideration per- 
formed not on the factor level but on the level of correla- 
tion functions This points to the fact that the result 
does not depend on the value of the density parameter 
of the system. 

Thus, the surface tension corresponds exactly to the 
number density oscillations (deviations) near the surface. 
We emphasize again that we are still talking about the 
boundary of the fluid and hard solid. 
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